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I.  Introduction 

The  work  for  this  research  contract  on  atmospheric-pressure 
gas  lasers  had  as  its  main  objectives  the  generation  of  short 
laser  pulses  through  mode  locking  and  cavity  dumping,  the  study 
of  the  nonlinear  amplification  processes  in  a transversely  excited 
atmospheric  pressure  laser,  and  development  of  computer  programs 
to  check  the  observed  amplifier  responses  against  models  of  the 
relaxation  processes  of  the  C02  molecule  in  a laser  discharge. 

Thus,  the  proposed  investigation  concentrated  on  three  methods 

by  which  pulses  could  be  made  both  more  energetic  and  of  shorter 
duration : 

(a)  mode  locking 

(b)  cavity  dumping 

(c)  nonlinear  amplification . 

The  work  on  nonlinear  processes  in  CC>2  systems  was  origin- 
ally begun  in  order  to  demonstrate  the  possibility  of  amplifying 
pulses  of  duration  short  compared  with  the  inverse  bandwidth  of 
the  amplifying  medium,  without  the  undesirable  pulse-broadening 
associated  with  linear  amplification.  Obviously,  this  can  be  done 
only  if  the  effective  bandwidth  of  the  nonlinear  amplification  is 
broader  than  that  of  linear  amplification.  Line  broadening  due 
to  optical  saturation  by  cw  radiation  is  well  known.115  The 
purpose  was  to  demonstrate -that,  dynamically,  line  broadening  via 
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saturation  could  be  utilized  to  amplify  over  a broader  bandwidth 
than  that  accorded  to  linear  amplification.  It  was  also  hoped 
that  dynamic  pulse  shaping  could  be  utilized  to  lead  to  additional 
pulse  shortening.  This  dynamic  pulse  shaping  is  closely  related 
to  the  optical  "nutation " effect.  * 3 > The  principle  of  pulse  shaping 
and  shortening  by  optical  nutation  was  demonstrated f4 1 just  before 
the  initiation  of  the  present  contract. 

Hie  work  under  this  contract  was  concerned  with  the  improve- 
ment of  the  pulse  shaping  and  shortening,  and  with  the  effort  to 
increase  energy  extraction  in  the  nonlinear  amplification.  The 
mode  locking-cavity  dumping  scheme  to  produce  single,  high  power 
pulses  was  abandoned  in  favor  of  a pulse  selection  outside  the 
cavity.  It  turned  out  that  the  ultimate  power  achievable  by  the 
pulse  selection  outside  the  cavity  Was  higher  than  the  one  achieved 
through  cavity,  doping.  The  electrooptic  crystal  used  for  cavity 
dumping  depressed  the  peak  power  of  the  laser  more  than  the  loss 
of  power  due  to  the  partial  mirror  transmission,  with  a conven- 
tional TEA  laser,  pulses  of  500  kw  peak  power  and  2 ns  long 
were  obtained.  These  were  used  as  the  experimental  probe  for 
short  pulse  amplification  studies. 

These  studies  were  aimed  at  determining  the  limits  on  energy 
extraction  and  achievable  pulse-widths  in  the  pulse  generation 
and  nonlinear  amplification.  One  problem  encountered  in  the 
amplification  of  short  pulses  in  COj  is  the  incomplete  energy 
extraction  from  the  rotational  levels  if  the  pulse  duration  becomes 
comparable  to  the  rotational  relaxation  time. 
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In  fact  the  extraction  time  of  the  energy  from  all  rotational 
* 5 ] 

levels  is  roughly  Neff  Trot>  where  Trot  the  ro-ational 

relaxation  time  and  N cc  is  the  effective  number  of  rotational 

err 

levels  (a  function  of  temperature,  a number  of  the  order  of  10)  . 
One  of  our  studies  concentrated  on  the  utilization  of  the  energy 
stored  in  the  vibrational  levels  as  discussed  in  Section  I. 

In  the  analysis  of  the  data  of  Section  I,  we  used  a computer 
program  that  treated  the  rotational-vibrational  relaxation  of  C02 
under  short-pulse  excitation.  This  computer  program  was  developed 
by  A.  Ross  of  this  laboratory  who  was  supported  by  a NSF  contract. 
The  computer  program  used  curves  obtained  by  Nighan^  for  the 
pumping  rates  of  the  energy  levels  by  the  discharge.  Nighan's 
work  is  based  on  a full  treatment  of  the  Boltzmann  equation  and 
the  evaluation  of  the  pumping  rates  requires  an  elaborate  computer 
program  and  considerable  computer  time.  We  considered  it  expedient 
to  develop  an  analytic  theory  of  electronic  pumping  rates  using 
approximate  cross  sections  for  the  vibrational  and  electronic 
excitation  rates.  The  advantage  of  this  approach  is  that  it  lends 
insight  into  the  important  mechanisms  determining  laser  efficiency. 
A description  of  the  theory  forms  section  II,  details  appear  in 
the  literature.^ 

The  study  of  energy  extraction  by  short  pulses  showed  the 
importance  of  obtaining  short  and  energetic  pulses  by  the  mode- 
locking process  in  the  first  place.  Saturable  absorber  modelockjLng 
has,  thus  far,  produced  much  shorter  mode  locked  pulses  than  forced 
modelocking  in  those  systoms  in  which  both  methods  have  been 
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employed  successfully.  Therefore  we  looked  seriously  at  node- 
locking  of  the  C02  system  by  saturable  absorbers,  both  experi- 
mentally and  theoretically.  We  found  that  no  adequate  analytic 
theory  of  saturable  absorber  mode  locking  has  been  published  in 
the  literature.  Therefore,  we  set  out  to  develop  such  a theory 
with  the  aim  of  attaining  a description  of  saturable  absorber 
modelocking  as  complete  as  that  presented  by  Kuizenga  and  Siegman[8J 
for  forced  modelocking.  We  were  successful  in  doing  this.  We 
believe  that  the  theoretical  work  on  saturable  absorber  modelocking 
is  the  most  important  success  of  the  work  performed  under  this 

contract.  Three  papers  on  the  subject  of  modelocking  have  appeared, 
or  will  appear  shortly.  15.10,111  The  1#8t  of  was  ^ 

at  Bell  Laboratories  by  the  principal  investigator,  who  spent  a 
year  with  the  group  working  on  modelockin  j of  dye  lasers  (Prof. 

Szflke  took  his -place  during  his  absence) . a further  paper  has 
been  submitted  for  publication.  It  forms  Appendix  I.  We  shall 
not  reproduce  the  papers  here  but  give,  in  Section  III,  an  account 
of  the  salient  features  of  the  theory  and  its  success  in  explaining 
certain  experimental  observations  on  passively  modelocked  systems. 

The  experimental  work  on  saturable  absorber  modelocking  con- 
centrated on  the  C02  system  using  SFg  as  a saturable  absorber. 

In  order  to  describe  the  action  of  the  SFg  absorber  adequately, 
a model  of  SF6  was  developed  which  accounts  better  for  the  trans- 
ient absorption  of  SFg  as  a function  of  pulse  energy  and  pulse- 

width  than  previously  proposed  models  . 112 '13J  This  topic  is  taken 
ap  in  Section  IV. 
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L: Energy  Extraction  from  Non-Lasing  Vibrational  Levels 

As  mentioned  in  the  introduction,  we  found  it  expedient  to 
select  modelocked  pulses  outside  the  cavity,  rather  than  via 
cavity  dumping,  because  the  final  energy  of  the  pulse  was  higher 
with  this  method  of  pulse  selection.  Also,  in  this  way  more  than 
one  pulse  could  be  "clipped"  from  the  pulse  train.  Two  pulses  in 
succession  were  used  to  probe  the  vibrational  population  recovery 
in  a C02  TEA  amplifier  and  hence  the  possibility  of  additional 
energy  extraction  from  the  amplifier,  within  times  corresponding 
to  the  vibrational  relaxation  times. 

The  amplifier  was  filled  with  a mixture  of  CC>2 , N2  and 
helium  at  200  Torr  and  the  partial  pressure  of  CC>2  was  varied 
in  the  experiment  from  5 Torr  to  40  Torr.  The  individual  input 
pulse  lengths  were  2 ns,  the  spacing  between  the  pulses,  12  ns. 
The  first  laser  pulse  was  made  to  enter  the  amplifier  ^30  ys 
after  the  application  of  the  discharge  current  pulse.  The  re- 
covery of  the  gain  of  the  second  pulse  as  a function  of  partial 
pressure  of  C02  was  determined.  The  total  pressure  in  the 
amplifier  was  sufficient  for  appreciable  equilibration  of  the 
rotational  population  distribution  within  the  duration  of  one 
pulse.  Furthermore,  since  the  rotational  relaxation  is  mainly  a 
function  of  the  total  pressure,  and  not  of  the  partial  pressure 
variation  of  the  partial  pressure  did  not  affect  the 


of  C02, 
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recovery  time  of  the  rotational  population  distribution.  There- 
fore the  increased  recovery  as  a function  of  increasing  pressure 
is  attributable  to  V-V  relaxation  processes.  The  level  which 
has  an  appreciable  population  after  equilibration  of  the  vibra- 
tional temperatures  and  which  can  feed  the  upper  laser  level  ap- 


preciably within  12  ns  at  the  partial  pressure  used  is  the  Oil 
level.  Hence  this  experiment  ascertains  the  rate  of  relaxation 
of  the  Oil  level  and  the  degree  of  population  recovery  by  means 
of  this  relaxation  process.  Figure  1 shows  the  result  of  the 
experiment.  At  the  very  low  partial  pressures,  when  the  relaxation 
mechanism  is  not  operative  within  12  ns,  the  gain  of  the  second 


pulse  compared  with  that  of  the  first  one  has  been  reduced  by  20% 
because  of  the  population  depletion  of  the  first  pulse.  The  ex- 
perimentally observed  recovery  as  a function  of  partial  pressure 
is  shown.  Two- theoretical  curves  computed  on  the  assumption  that 
th-  recovery  is  due  entirely  to  the  feeding  of  the  upper  laser 
level  by  the  Oil  level  population  are  shown  as  dashed  curves. 

The  two  curves  were  computed  for  two  values  of  population  inversion 
changes  consistent  with  the  gain  changes  caused  by  the  first  pulse. 
This  shows  that  the  major  portion  of  the  recovery  is  explainable 
by  this  mechanism.  The  fact  that  the  experimentally  observed 
recovery  is  larger  than  the  one  predicted  on  the  basis  of  the  Oil 
level  relaxation  is  attributable  to  the  more  effective  relaxation 
of  the  lower  laser  level  with  increasing  partial  pressure,  an 
effect  not  taken  into  account  in  the  theory  which  assumed  that 


even  at  the  low  partial  pressures  the  lower  laser  level 


relaxes 
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completely  within  the  12-ns  interval. 

Details  of  the  experiment  are  given  in  the  RLE  Quarterly 
Progress  Report  No.  Ill,  Oct.  15,  1973.  ».  theoretical  compu- 

tations are  contained  in  the  Master's  thesis  of  V.  Manichaikul . 
The  conclusion  drawn  from  the  experiment  is  that  energy  storages 
available  in  vibrational  combination  modes  can  be  utilised  if 
the  pulse  to  be  amplified  is  lengthened.  The  present  experiment 

indicates  the  magnitude  of  the  effect  and  the  time  scales  at  which 
such  utilization  occurs. 
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Electron  Distribution  and  Lasing  Effiripn^v 


Computer  codes  evaluating  the  electron  distribution  and  the 
pumping  of  laser  levels  have  been  developed  in  several  laboratories 
notably  at  United  Aircraft  Corporation . [15 1 Whereas  good  quantita- 
tive results  can  be  obtained  from  these  programs,  it  is  difficult 
to  cull  from  them  physical  insight  about  the  influence  of  discharge 
parameters  on  gain  and  lasing  efficiency,  without  carrying  out 
costly  computations.  With  this  in  mind,  we  have  developed  a sim- 
plified model  for  the  electron  distribution  and  the  pumping  of  the 
electron  laser  levels  which  is  amenable  to  closed-form  solutions. 
The  details  of  this  work  have  been  published.171  Here  we  give  a 
summary  of  the  salient  features. 

The  energy  transfer  of  the  electrons  to  the  vibrational  levels 
and  to  the  electronic  radiation  levels  (which  are  useless  energy 
sinks  in  the  case  of  a IR  molecular  laser)  is  controlled  by  the 
shape  of  the  velocity  distribution  function  of  the  discharge 
electrons.  This  distribution  function  is  approximately  spherically 
symmetric  and  hence  describable  in  terms  of  a function  f (v)  of 
a scalar  v,  the  magnitude  of  the  electron  velocity. 

The  basic  equation  governing  fo(v)  is  derived  from  the 
Boltzmann  Equation: 


where  vc  is  the  frequency  of  momentum  transfer  collisions  and 
E is  the  electric  field  in  the  discharge.  The  term  on  the  left- 
hand  side  represents  the  rate  of  entry  of  electrons  into  the 
velocity  range  v,  v + dv  produced  by  the  driving  electric  field. 
The  rate  of  entry  of  electrons  is  equal  to  the  rate  of  exit  caused 
by  the  inelastic  collisions,  the  term  on  the  right-hand  side  of 
the  equation.  The  distribution  fuction  fQ(v)  is  normalized: 

J 4ttv2  fQ(v)  dv  = 1. 

An  analytic  treatment  results  from  an  appropriate  modeling  of  the 
inelastic  collision  term  on  the  right  hand  side  of  (1) . One  ap- 
proximates the  collision  cross  sections  of  the  electrons  with  the 

molecules,  which  exhibit  resonant  peaks  as  functions  of  energy, 

* 

by  impulse  functions  of  energy.  The  main  collision  processes  are 
the  vibrational  energy  excitations  with  peaks  near  2 eV  (for 
N2  and  C02),  and  electronic  excitations,  with  peaks  that  may  all 
be  lumped  into  a single  excitation  cross  section,  as  long  as  the 
average  electron  energies  are  as  low  as  those  generally  produced 
in  molecular-laser  gas  discharges. 

With  these  assumptions,  closed  form  solutions  can  be  ob- 
tained for  fQ ( v)  (or  fQ(u)  where  u = mv2/2c)  . 

Some  interesting  conclusions  were  obtained  from  the  analysis, 
(a)  In  an  E-beam  excited  discharge  we  find  that  the  V-I 
characteristic  of  the  sustainer  electrodes  remains  unchanged  by 
the  presence,  or  absence,  of  lasing.  For  a given  power  input 
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per  unit  volume,  lasing  redistributes  the  flow  of  energy  among 
the  levels,  changes  the  electron  distribution,  but  does  not  affect 
the  conductivity  of  the  medium.  This  conclusion  turned  out  to  be 
rather  model  independent.  The  sole  source  for  this  somewhat  sur- 
prising result  is  the  assumption  that  the  elastic  collision 
frequency  is  energy  independent. 

(b)  The  efficiency  of  a E-beam  laser  peaks  at  a value  of 
electric  field  that  maintains  the  average  energy  of  the  electrons 
well  below  2 eV,  the  energy  at  which  the  vibrational  excitation 
cress  section  was  assumed  to  peak.  The  reason  for  the  rapid  de- 
crease of  efficiency  with  increasing  electric  field  beyond  this 
critical  field  is  that  the  energetic  "tail”  of  the  electron  dis- 
tribution excites  the  electronic  levels  excessively. 

While  these  conclusions  could  have  been  culled  from  Nighan's 
* 

computer  analysis,  it  is  helpful  that  the  analytic  treatment  reaches 
these  conclusions  in  a transparent  way  and,  in  fact,  gives  results 
even  in  (approximate)  quantitative  agreement  with  Nighan's  results. 
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III.  Theory  of  Passive  Modelocking  of  a Homogeneously 

Broadened  Laser 

The  purpose  of  the  new  theory  of  passive  modelocking 10 
is  to  provide  an  analytic  description  of  the  phenomenon  equivalent  to 
the  one  presented  by  Siegman  and  Kuizenga^  for  the  case  of  active 
modelocking.  Passive  modelocking  is  a phenomenon  more  complicated 
than  active  (forced)  modelocking  because  the  modulation  provided 
by  the  saturable  absorber  is  a function  of  the  pulse  amplitude 
and  shape  passing  through  the  absorber,  and  depends  on  the  relaxa- 
tion time  of  the  absorber. 

Just  as  in  the  case  of  active  modelocking,  the  theory  ot 
passive  modelocking  is  concerned  at  first  with  cw  operation,  because 
it  is  simpler.  We  shall  first  outline  the  theory  and  then  present 
arguments  that  many  results  of  the  theory  are  applicable  to  the 
transient  buildup  of  modelocking. 

The  analysis  makes  the  following  assumptions  and  approximations: 

(a)  The  laser  medium  is  homogeneously  broadened. 

(b)  A pulse  passing  through  any  one  of  the  elements  of  the  system, 
saturable  absorber,  gain  medium,  etc.,  is  modified  only 
slightly  (say  20  percent  gain  or  loss)  on  one  single  pass. 

This  assumption  enables  one  to  expand  exponentials  to  first 
order  in  their  argument. 

(c)  The  dispersion  of  the  system  as  a function  of  frequency  may 
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be  expanded  to  second  order  in  frequency  (this  assumption 
is  legitimate  if  the  laser  pulse  spectrum  is  narrow  compared 
with  the  bandwidth  of  the  system)  . 


Under  these  assumptions  one  arrives  at  the  following  equation  for 
the  envelope  v(t)  of  the  electric  field,  E(t)  = v(t)  exp  ju  t, 

of  the  pulse  traveling  in  one  direction  between  the  mirrors  of  the 
the  laser  cavity: 


♦ jb  - g 


+ 6 + q d 


u>. 


dt 


m(t) 


v(t)  = 0 


(1) 


Here  the  factor  1 expresses  the  effect  of  the  normalized  cavity 
loss  on  the  pulse,  jb  is  the  reactive  effect  caused  by  oscillation 
of  the  cavity  modes  off  their  natural  resonance  frequency; 

is  the  operator  describing  the  effect  of  the 

dispersive  laser  medium  gain  on  the  pulse  envelope;  a>L  is  the 
laser  medium  linewidth  and  g is  the  gain  at  line  center  per  pass 
normalized  to  the  loss  per  pass.  The  operator  d2/dt2  is  a dif- 
fusion operator  in  "time-space"  which  accounts  for  the  spreading 
of  the  pulse  in  time  because  of  medium  dispersion.  The  operator 

~ 2 IT  Ascribes  net  tirie  advance  5T  (or  delay  -6T  of 
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the  pulse)  with  respect  to  the  empty  cavity  roundtrip  time  T 

R 

of  the  cavity  as  caused  by  the  laser  medium  dispersion,  as  well 
as  that  induced  by  the  modelocking  modulation  m(t)  of  the  node- 
locking  element.  Equation  (1)  sets  equal  to  zero  the  sum  of  the 
modifications  of  the  pulse  due  to  cavity  loss,  laser  gain,  pulse- 
advance  and  modulation,  i.e.  it  states  the  requirement  of  a steady 
state.  The  fundamental  relation  (1)  describes  both  the  forced, 

and  passive  methods  of  modelocking.  Indeed,  if  the  nodelocking  is 
forced,  then 

m(t)  = 2M [1  - cos  wmt] 

where  M expresses  the  modulation  amplitude  normalized  to  the 
cavity  loss,  and  is  the  frequency  of  modulation.  When  (2) 

is  introduced  into  (1)  and  after  the  change  of  variable  v(t)  = V(t) 
exp  (6  + g/2)  ioLt,  one  obtains  an  equation  that  is  identical  in 
form  to  the  equation  of  notion  of  quantum  mechanics  of  a particle 
in  a periodic  potential  (the  time  t must  be  replaced  by  the  spa- 
tial coordinate  x) . The  result  of  Siegman  and  Kuizenga  is  ob- 
tained through  one  further  approximation.  If  the  modulation  is 
strong,  leading  to  well  isolated  pulses,  then  1 - cos  co  t may  be 
approximated  by  1/2  (wmt) 2 within  the  time  interval  occupied  by 
the  pulse.  Under  this  approximation  the  well  known  gaussian  pulse 
solution  is  obtained,  in  addition  to  the  gaussian,  higher  order 
Hermite  gaussian  solutions  are  obtained  as  well  as  previously 
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pointed  out  by  Haken  and  Pauthier. [14J  We  were  able  to  show[9] 
that  all  pulse  solutions,  except  the  lowest  order  Hermite  Gaussian, 
are  unstable  for  synchronous  modulation  («m  is  set  equal  to  the 
2ir/Tp  where  Tp  is  the  cavity  roundtrip  time  including  the  ad- 
ditional delay  caused  by  dispersion  of  the  laser  medium) . 

Hence,  Equation  (1)  describes  fully  the  forced-modelocking 
results  of  Siegman  and  Kuizenga.  Saturable  absorber  modelocking 
is  described,  when  the  modulation  m(t)  is  determined  self  con- 
sistently with  the  pulse  intensity  passing  through  the  saturable 

absorber.  For  an  absorber  with  relaxation  time  fast  compared  to 
the  pulse  duration, 

m(t)  « i a i - H2 

(■ : £)  '' 

where  PA  expresses  the  saturation  power  of  the  absorber,  if  v 
is  so  normalized  that  | vj 2 is  equal  to  the  instantaneous  power 
in  the  pulse.  The  resulting  nonlinear  differential  equation  can 
be  solved.  The  solutions  are  elliptic  integrals  which  reduce  to 
a secant  hyperbolic  in  the  limit  of  well  separated  pulses.  [10] 

If  the  absorber  is  slow  compared  with  the  pulse  duration,  then 
one  faces  a stability  problem:  if  the  absorber  does  not  "shut-off" 

after  passage  of  the  pulse,  perturbations  following  the  pulse  will 
grow  and  the  solution  is  unstable.  Stable  pulses  can  be  obtained 
xf  the  laser  gain  is  reduced  by  the  passage  of  the  pulse  so  as  to 
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shut  off  the  net  gain  after  the  pulse  has  passed.  Again,  under 
certain,  not  too  restrictive  assumptions,  the  pulse  solutions  are 
secant  hyperbolics . It  can  be  shown  in  general  that  the  wings 

of  pulses  produced  by  saturable  absorber  modelocking  must  be 
exponential,  not  gaussian . This  is  of  importance  when  using 
_hese  pulses  in  nonlinear  time  resolved  spectroscopy  on  a time 
scale  comparable  to  the  pulse  duration.  This  exponential  time 
dependence  of  the  pulse  wings  has  been  observed  experimentally  on 
modelocked  dye  laser  pulses. 

So  far  we  have  discussed  theories  of  cw  passive  modelocking 
We  contend  that  the  theories  are  relevant  to  the  transient  build- 
up of  modelocking  provided  certain  criteria  are  met.  We  further 
claim  that  these  criteria  must  be  met  if  the  transient  modelocking 
of  a gain  switched  laser  is  to  provide  deterministic  (nonstatistical) 
pulse  trains,  repeatable  from  "shot"  to  "shot".  We  shall  come  back 
to  this  point  after  a brief  discussion  of  previous  work  on  the 
transient  buildup  of  modelocked  pulses. 

Much  computational  work  was  done  in  the  literature  ^ 

the  transient  buildup  of  saturable  absorber  modelocking.  The  pre- 
vailing view  was  the  one  originally  proposed  by  Letokhov : 

The  buildup  of  modelocked  pulses  is  started  by  initial  noise  spikes 
of  random  width  and  height  traveling  between  the  two  laser  reflec- 
tors. From  among  these  spikes,  those  of  largest  energy  are  selected 
by  the  "gating  action"  of  the  saturable  absorber.  A natural  con- 
sequence of  this  model  is  that  the  predicted  pulsed  output  of  a 
Q-switched  or  gain  switched  laser  is  statistical  in  nature. 
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i , 

Experimen tally , no  doubt/  the  modelocked  pulse  from  gain 
switched  rJd  Glass  or  Nd  YAG  do  exhibit  statistical  behavior. 

This  can  be  undesirable  if  predictable  pulses  are  to  be  produced 
from  shot  to  shot.  Special  experimental  efforts  have  been  made 
to  produce  repeatable  pulse  trains  of  modolockod  pulses.  This 
has  always  been  found  to  necessitate  operation  of  the  laser  near 
threshold.  It  is  our  contention  that  there  is  a limit  of  transient 
modelocked  operation  in  which  the  modeiocking  is  entirely  determin- 
istic. The  steady  state  theory  of  modelocking  bears  on  the  build- 
up of  these  "deterministically"  modelocked  pulses. 

In  support  of  this  picture  we  cite  the  following  argument. 
Initially,  the  laser  builds  up  from  spontaneous  emission  no^'se. 

In  order  to  reach  power  levels  that  cause  saturation  of  the  absorber, 
the  power  has  to  grow  by  a factor  of  1015  or  so.  Such  large 
growth  leads  to  gain-bandwidth  narrowing  that  is  proportional  to 
the  logarithm  of  the  growth  factor — an  appreciable  factor.  Hence 
the  power  traveling  in  one  direction  inside  the  cavity  is  smoothed 
by  the  bandwidth  narrowing.  All  of  this  still  occurs  during  the 
buildup  process.  Of  all  Fourier  components — harmonics  of  2m/T 

R 

the  lowest  Fourier  component  grows  fastest  initially.  Hence 
the  'smoothed"  average  (with  some  statistical  structure)  will 
evolve  into  a growing  sinusoidal  perturbation  of  a "dc  background". 

If  the  system  is  not  driven  excessively  above  threshold,  then  the 
time  evolution  of  the  pulse  may  be  treated  adiabatically , as  a 

succession  of  quasisteady  states.  This  theory  is  outlined  in 
Appendix  II. 
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— — Short  Pulse  Saturation  of  High  Pressure  SF_ 

— ■ — ——6 

Interest  in  the  short  pulse  saturation  properties  of  SF 
has  been  stimulated  by  its  usefulness  as  a Q switching [26 1 and 
modelocking [27]  element  for  the  P branch  of  the  10.4  p band 
of  the  C02  laser.  More  recently,  work  on  laser  isotope  separa- 
tion using  SFg  has  intensified  that  interest/281 

Our  primary  motivation  in  carrying  out  the  present  study  of 
SF6  short  pulse  saturation  was  to  determine  the  characteristics 
of  the  medium  appropriate  to  nodelocking  the  high  pressure  a>2 
laser.  The  modelocking  mechanism  is  the  bleaching  of  the  SFg 
absorption  on  the  time  scale  of  the  modelocked  pulses.  Thus,  the 
absorption  is,  in  part,  due  to  a set  of  levels  having  a relaxation 
time  on  the  order  of  or  less  than  the  typical  width  of  the  mode- 
locked  pulse.  . In  other  words,  the  attenuation  experienced  by  an 
intense  pulse  of  duration  Tp  in  propagating  through  a cell  of 
SF6  should  be  less  than  that  experienced  by  a quasi-continuous 
signal  of  the  same  energy.  By  quasi-continuous  we  mean  a pulse 
of  duration  Tp  greater  than  that  of  the  modelocked  pulse  but 

less  than  the  V-T  relaxation  time  of  SF  . 

6 

We  have  completed  a set  of  experiments  which  verify  this 
hypothesis.  By  alternately  irradiating  a 1 mm  long  cell  of  SF 

6 

with  the  single  mode  output  and  the  modelocked  output  of  a high 
pressure  TEA  CC>2  laser  we  have  measured  the  difference  in  trans- 
mission of  the  cell  for  equal  values  of  the  total  energy  in  the 
TEA  pulse.  Single  mode  operation  was  achieved 
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via  a 2.5  cm  long  intracavity  cell  of  low  pressure  SF  (<  1 
torr)  . Modelocking  was  produced  by  an  intracavity  acousto-optic 
modulator  tuned  to  the  cavity  mode  spacing  at  approximately  42 
MHz.  The  mode locked  pulses  were  2 ns  long  spaced  24  ns  apart, 
thus  giving  approximately  a 12:1  peak  power  enhancement  over 
the  single  mode  TEA  pulse  of  equal  energy. 

The  output  of  the  TEA  laser  was  passed  through  a 3 mm 
aperture  to  ensure  a nearly  uniform  intensity  distribution  across 
the  beam  and  focussed  by  a salt  lens  to  a .2  mm  diameter  spot 
inside  the  SFg  cell,  a second  lens  recollimated  the  beam  before 
it  was  detected  by  a Au:Ge  detector.  Different  values  of  input 
intensity  were  achieved  by  placing  calibrated  CaF2  attenuators 
in  the  path  of  the  beam.  From  the  detector  the  signal  was  ampli- 
fied by  an  amplifier  which  was  slow  (-v  5ys)  compared  to  the 
TEA  pulse  ('V/  '250  ns)  so  that  its  response  depended  only  on  the 
total  energy  in  the  TEA  pulse.  The  amplifier  output  was  then 
fed  into  a RIDL  400  channel  pulse  height  analyzer  which  displayed 
the  distribution  of  transmitted  TEA  pulse  energy  for  successive 
shots  .of  the  laser.  The  use  of  the  analyzer  enabled  us  to  dis- 
tinguish between  shot  to  shot  fluctuations  in  the  laser  output 
energy  ('v,  10%)  and  differences  in  absorber  transmission. 

Results  of  the  transmission  measurements  at  various  SFr 
pressures  for  the  P(22)  line  of  10.6  y C02  output  are  shown 
in  Figure  2.  The  modelocked  pulse  transmission  is  always  higher 
than  that  of  the  single  mode  pulse.  Furthermore,  we  observe  that 
the  saturation  curves  appear  to  be  roughly  asymptotic  to  the  small 
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signal  absorption  values  measured  by  Wood  et  al.^29^  indicating 
that  the  absorption  cross-section  does  not  change  appreciably  with 


pressure . 

The  complexity  of  the  energy  level  structure  of  SF  makes 

6 

it  difficult  to  analyze.  Models  developed  by  Burak  et  al.fl2J  and 
Brunet  show  excellent  agreement  with  experiments  on  the  cw 
saturation  of  SF6 . The  response  of  SF,.  to  short  pulses,  how- 
ever, is  not  well  understood.  Recent  work  by  Oppenheim  and 
Kaufman  shows  that  the  pulse  saturation  behavior  of  SF  can 
be  dramatically  different  from  the  cw  behavior.  However,  the 
model  adopted  by  Oppenheim  and  Kaufman  is  unrealistic  in  that  it 
assumes  an  instantaneous  relaxation  rate  for  the  excited  state 
transition.  Pulse  transmission  measurements  have  also  been  per- 
formed by  Armstrong  and  Gaddy,  !31)  but  their  model  ignores  V-V 
equilibration  and  neglects  the  v3  - 2v3  excited  state  absorption. 
The  model's  prediction  of  a vg  - Vj  + Vf.  relaxation  time  which 
is  much  faster  than  gas  kinetic  and  a V-T  time  which  is  three 
orders  of  magnitude  faster  than  previously  determined  values'321 


.is  a measure  of  its  weakness. 

We  have  developed  a simple  closed  form  theory  to  explain  the 
saturation  curves  for  both  the  single  mode  and  modelockod  pulses. 
Our  approach  is  similar  to  that  of  Oppenheim  and  Kaufman  in  that 
we  consider  only  the  absorption  of  the  v3  mode  and  neglect  the 
coupling  to  other  vibrational  modes  of  the  molecule.  In  contrast 
to  previous  work,  however,  we  include  the  three  fold  degeneracy 
and  full  V-V  equilibration  of  the  Vj  mode  and  the  saturation 
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of  the  excited  state  absorption  v3  - 2v3.  To  do  so  we  adopt  the 
model  shown  in  Figure  3 in  which  the  following  notation  is  used: 


Nl'  N2'  N2  ’ ' N3  Population  densities  in  the  three  vibrational 

levels  which  are  directly  interacting  with  the 
incident  field. 


Nir  N22'  N33  Population  densities  which  do  not  interact 

directly  with  the  field. 

V °l  Absorption  cross-sections  of  the  ground  state 

and  first  excited  state  respectively. 


Rotational  relaxation  rate. 

Bv  Fraction  of  total  number  of  molecules  in  the 

vth  vibrational  level  which  interact  directly 
with  the  field. 


k 


w 


Rate  constant  for  equilibration  among  the 
vibrational  levels. 


We  have  omitted  any  reference  to  the  V-T  relaxation  time  because 
it  is  slow  compared  to  our  pulse  duration.  Measured  values  of  the 
V-T  rate  place  it  at  122  ps/torr1321  which  is  much  greater  than 
our  250  ns  pulsewidth  even  at  250  torr,  the  maximum  pressure 
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used  in  our  experiments. 

The  theoretical  saturation  curves  computed  from  the  above 
model  are  compared  to  the  experimental  single  mode  P(22)  trans- 
mission data  in  Figure  4 . Good  agreement  is  found  at  low  pressures 
At  higher  pressures,  however,  SFg  is  harder  to  saturate  than  the 
model  predicts.  This  discrepancy  can  be  attributed  to  the  coupling 
between  and  other  vibrational  modes  of  the  i olecule  which  is 

neglected  in  our  model.  The  difference  between  theory  and  exper- 
iment indicates  that  approximately  three  vibrational  modes  con- 
tribute to  the  absorption  at  pressures  of  100  torr  and  greater. 

The  details  of  this  work  will  be  submitted  for  publication. 
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Conclusions 

Of  the  results  obtained  by  the  research  summarized  here 
the  most  important  one  is  believed  to  be  the  development  of  a 
theory  of  saturable  absorber  modelocking.  In  the  process  of 
its  application  to  C02  laser  modelocking  a better  understanding 
of  SF6  as  a modelocking  absorber  has  been  obtained.  Numerous 
publications  have  resulted  from  the  work,  some  bearing  specific 
acknowledgement  to  this  contract,  some  which  have  been  conceived 
as  an  outgrowth  of  this  work  while  the  principal  investigator 
was  on  a leave  of  absence  at  the  Bell  Laboratories,  in  part  with 
the  objective  of  applying  the  new  theory  of  saturable  absorber 
modelocking  to  the  cw  modelocked  dye  laser  producing  subpico- 
second pulses  at  the  Bell  Laboratories,  a list  of  publications, 

which  grew  out  of  this  work,  published,  accepted,  and/or  submitted 
is  summarized  below: 

With  Contract  Acknowledgment  Published: 

li]  W.  p.  Allis,  H.  A.  Haus,  "Electron  Distributions  in  Gas 
Lasers,"  J,  Appl . Phys . , £5,  781-791  (1974). 

lii]  H,  A,  Haus,  "A  Theory  of  Forced  Modelocking,"  IEEE  J. 

Quantum  Electronics.  QE-11,  323-330  (1975) 

liii]  H.  A.  Haus,  "Theory  of  Modelocking  with  a Fast  Saturable 
Absorber,"  J.  Appl.  Phys..  40 , 3049-3059  (1975) 

" 
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Submi  tted: 

Uv]  p.  Hagelstein,  C.  P.  Ausschnitt,  "Shape  and  Stability 
Dependence  of  Passively  Model ocked  Pulses  on  Absorber 
Relaxation  Time,"  to  J.  Appl.  Phvs . 

Without  specific  contract  acknowledgment,  because  carried  out 
at  Bell,  but  conceived  as  an  outgrowth  of  the  present  work: 

To  be  published: 

[v]  H.  A.  Haus,  "Theory  of  Modelocking  with  Slow  Saturable 
Absorber,"  to  be  published  in  J.  Quantum  Electronics. 
September  (1975) 


(Vi)  H.  A.  Haus,  C.  v.  Shank,  E.  P.  Ippen,  “Shape  of  Passively 
Mode locked  Laser  Pulses,"  Optics  Comm. 


[viij  II.  a.  Haus,  "Model  of 
Appl . Phys . 


cw  Modelocked  Dye  Laser,"  to  J. 


[viii]  H.  A.  Haus,  Parameter  Ranges  for  cw  Passive  Mode  lock  i n g,  " 
to  IEEE  J.  Quantum  Electronics . 
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Personnel  Working  on 


Students : 


Yongyut  «anichi*ul  (Master,s  Thesis  ^ 

and  Amplification  of  High  Intensity  Nanosecond  raA  Laser 

Pulses  , " September  1973), 


Peter  Hagelstein  (Bachelor’s  Degree,  Thesis  title:  ’.Mode locking 

by  saturable  Absorb,  r with  ^negligible  Aelaxation  Tims,’,  to 
be  submitted  in  September  19  75) 


Christopher  Ausschnitt  (Doctor’s  Degree,  Thesis  title:  ’’Hode- 

locking  Studies  on  High  Pressure  CO,  haters ,’’  to  be  submitted 
in  September  1975) . 


Faculty; 


H.  A.  Haus 


A.  Szflke  (took  Professor  Haus'  place  duri 
Bell  Laboratories) 


ng  leave  of  absence  at 


W.  Allis* 


Not  supported  financially. 
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Abstract  ? 

The  dependence  of  the  saturable  absorber  modelocked  pulse 
shape  upon  the  relaxation  time  of  the  absorber  is  determined  when 
the  laser  gain  is  assumed  constant.  The  pulses  are  found  to  be 
asymmetric  with  a slowly  rising  edge  and  more  rapidly  decaying 
edge.  The  modelocked  pulsewidth  is  a weak  function  of  absorber 
relaxation  time  TA,  becoming  broader  with  increasing  T . Fur- 
thermore, the  pulses  are  found  to  experience  a delay  in  traversing 
the  absorber  which  lengthens  the  cavity  transit  time.  The  growth 
of  perturbations  introduced  into  the  steady  state  is  shown  to  be 
bounded,  indicating  that  stable  modelocked  operation  is  possible 

even  in  cases  where  T^  is  comparable  to  or  greater  than  the  width 
of  the  steady  state  pulse. 
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Shape  and  Stability  Dependence  of  Passively  Modelocked 
Pulses  on  Absorber  Relaxation  Time 

P.  L.  Hagelstein , C.  P.  Ausschnitt 


Introduction 


Mode-locking  of  a laser  with  a passive  elementhas  been 
studied  by  many  authors.  Most  of  the  theoretical  studies  are 
either  computer  studies  ^ or  qualitative  analytic  studies  of 
pulse  shaping,  without  prediction  of  a pulse  shape  One 

exception  is  the  recent  work  by  Ilaus^3,4^,  who  obtains  closed  form 
expressions  for  the  pulse  shape  produced  by  passive  mode-locking 
in  two  limiting  cases:  (a)  The  relaxation  time  of  the  saturable 

absorber  is  short  compared  with  the  pulse  width,  (b)  the  relaxa- 
tion time  of  the  absorber  is  much  longer  than  the  pulsuwidth. 


In  case  (b)  the;  laser  gain  medium  must  contribute  to  the 
pulse  shaping^  '* ^ because,  otherwise,  the  system  would  possess 
gain  after  passage  of  the  pulse  and  perturbations  following  the 
pulse  would  experience  growth;  the  pulse  would  be  unstable. 

In  this  paper  we  study  the  intermediate  range  of  absorber 
relaxation  time,  when  the  absorber  relaxation  time  is  comparable 
to  the  pulsewidth . We  limit  our  analysis  to  the  case  where  the 
relaxation  time  of  the  gain  medium  greatly  exceeds  that  of  the 
absorber.  Since  there  may  be  net  gain  in  the  system  following 


r 


the  pulse  when  the  absorber  loss  is  relaxing  (in  the  case  of  well- 
separated  pulses  in  a system  with  long  absorber  relaxation  time) 
the  steady  state  pulse  found  may  be  unstable.  The  unstable  nature 
of  passive  mode-locking  with  a slow  absorber  in  the  case  where  the 
gain  is  hard  to  saturate  is  an  experimental  fact^^. 

Our  approach  is  similar  tc  that  taken  by  Creighton  et  alf7^, 

New  ^ , and  Haus  ^ in  that  it  focusses  on  modifications  of  a pulse 
in  one  round  trip  transit  through  the  laser  cavity.  The  steady  state 
is  obtained  by  requiring  the  pulse  to  reproduce  itself.  In  contrast 
to  Creighton,  the  saturable  absorber  modulation  is  included  self 
consistently.  The  gain  dispersion,  neglected  by  New  and  Letokhov 
in  the  nonlinear  regime,  is  taken  to  be  a major  factor  in  the  de- 
termination of  pulse  shape  and  width. 

In  section  I we  summarize  the  main  assumptions  and  derive 
» 

the  basic  equation.  Section  II  presents  some  features  of  the  pulse 
shape  that  can  be  gleaned  directly  from  the  basic  equation.  One 
can  show  that  the  pulse  shape  is  in  general  asymmetric  with  a grad- 
ually rising  edge.  In  Section  III  we  present  computer  solutions 
which  confirm  this  prediction  and  yield  the  two  main  pieces  of 
information  for  the  pulse  as  a function  of  Tfi,  its  width  and  its 
net  delay.  Finally,  in  Section  IV  we  investigate  the  stability 

of  the  solutions  and  find  that,  to  first  order,  perturbation  growth 
is  bounded. 


/ 
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I.  Assumptions  and  Basic  Equation 

The  system  considered  here  consists  of  three  elements  (see 
Figure  1) -the  laser  cavity,  the  active  gain  medium  and  the  satur- 
able absorber.  The  notation  and  method  of  derivation  parallels 
that  of  llaus[3].  The  following  assumptions  are  made: 

(1)  The  system  is  operating  in  the  steady-state  near  threshold. 

(2)  The  gain  medium  and  absorber  are  short  compared  to  the  cavity 
length  and  the  effect  of  each  element  on  a pulse  over  a single  pass 
is  small  (such  that  the  exponential  gain  and  loss  factors  can  be 
expanded  to  first  order)  . The  pulses  repeat  with  either  an  advance 
or  delay  from  the  cavity  transit  time,  the  magnitude  (assumed  small) 
and  sign  of  which  is  to  be  determined  by  requiring  the  pulse  to  re- 
produce itself  from  pass  to  pass. 

(3)  As  in  the  fast  absorber  mode  locking  solution  given  by  Haus, 
we  consider  only  the  case  where  the  pulse  repetition  time  is  much 
longer  than  the  pulsewidth . 

(4)  The  cavity  loss  is  constant  in  time  and  the  cavity  modes  are 
uniformly  spaced  in  frequency. 

(5)  The  active  medium  is  homogeneously  broadened  and  its  relaxa- 
tion time  is  much  longer  than  the  pulse  repetition  tine  (the  steady 
state  gain  is  approximately  constant  in  time) . 

(6)  The  lmewidth  of  the  gain  is  approximated  by  a Lorentzian  in 
frequency: 
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where  toQ  is  the  center  frequency  and  toL  is  a measur.  of  the 
bandwidth.  The  pulse  is  assumed  to  have  a bandwidth  small  compared 
to  the  lmewidth  so  that  the  gain  can  be  expanded  to  second  order 
in  frequency.  Transforming  into  the  time  domain,  the  frequency  de- 
pendence becomes  an  operator  by  replacing  (jw)n  by  (dn/dtn) 
wherever  it  occurs.  Thus,  the  gain  factor  operating  on  the  pulse 
envelope  v(t)  can  be  expressed  as  follows: 


exp 
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when  operating  on  the  pulse  envelope. 

(7)  The  bandwidth  of  the  absorber  is  wide  compared  to  the  ga 
linewidth . 


gam 


A two  level  quantum  system  with  saturation  and  relaxation  de- 
scribes the  absorber  population  density  difference. 


an 


(t)  _ n (t)  - ne  crA|  v(t)  i 2 n (t) 


dt 


•tito  A* 
o A 


(3) 


where  n(t)  is  the  papule  tier,  density  difference, 
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section  and  AA  is  the  beam  cross  section  for  the  absorber.  The 
degree  of  saturation  is  assumed  small  so  that  n(t)  = n in  the 
saturation  term  of  the  rate  equation.  The  relaxation  time  t 
is  comparable  to  a pulscwidth.  The  loss  due  to  the  absorber  is 

taken  to  be  L(t)  = aA  0A  AA  n (t)  where  0A  is  the  length  of  the 
absorber. 


A pulse  passing  through  the  reference  plane  (see  Figure  1)  on 

the  n-th  pass,  vn(t),  is  modified  after  making  a round  trip 

pass  through  the  cavity  such  that 
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where  exp  [- (uq/2Q) Tr]  is  the  constant  cavity  loss  factor  and 
exp [-2L ( t) ] is  the  time  dependent  absorber  loss  factor. 

t 

Requiring  the  pulse  to  be  periodic  in  time  and  expanding  the 

exponentials  to  first  order  results  in  the  following  constraint  on 
pulse  shape: 
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We  have  defined  the  same  quantities  normalized  to  the  cavity  loss 
used  by  Haus,  namely  the  small  signal  absorber  loss 


2Q 


the  saturated  laser  gain 
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the  absorber  saturation  power 
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and  the  time  shift  parameter 
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where  6T  is  the  advance  of  the  pulse  with  respect  to  the  empty 
cavity  round  trip  transit. 

The  first  term  on  the  left  of  (5)  is  the  normalized  constant  ca- 
vity loss.  The  next  two  terms  are  the  normalized  time -dependent  ab- 
sorber loss  where  q is  a measure  of  the  magnitude  of  the  loss.  The 
exponential  and  integral  are  the  solution  to  the  absorber  rate  equation 
(3)  . The  following  term  is  the  dispersive  gain  of  the  active  medium 
[(g/WjMd/dt)  acts  as  a delay  and  <g/Wj2  d2/dt2)  is  diffusive], 
where  g is  a measure  of  the  magnitude  of  the  saturated  gain.  The 
relative  magnitudes  of  the  unsaturated  losses  and  saturated  gain  are 
such  that  1 + q - g <<  g.  jn  other  words,  we  require  that  the  laser 
operate  near  threshold.  The  final  term  is  due  to  pulse  advance  or 
delay,  and  comes  about  from  a Taylor  expansion  in  time: 
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The  delay  or  advance 


must  be  determined  as  an  eigenvalue  of  (5) . 
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Equation  (5)  is  identical  to  (2.6)  of  Haus 


, except  that 


tho  response  of  the  saturable  absorber  is  now  represented  in  terms 

of  an  integral.  Tho  equation  can  bo  normalized  and  recast  in  a 
differential  form 
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9 dx  dx‘ 


y - y iyii  J2i 


dx  y dx2 


wt  T..  ^ 

- («  + q)  — y - 2y!  = 0 

g dx 


where 


y(x)  = 


v( t/r  ) 


- (1  + q - g) 
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The  time  is  normalized  to  tho  pulse  width  tp  of  the  fast  absorber 
™°del.  Tp  is  the  natural  time  scale  of  the  system.  For  a given  T , 
the  time  shift  parameter  ({  + g)  with  respect  to  the  cavity  round.*' 
trip  time  as  modified  by  the  laser  medium  is  determined  as  an  eigenvalue 
Therefore,  the  absorber  relaxation  time  T^  is  the  only  adjustable 
parameter  in  equation  (7).  In  the  limit  of  instantaneous  relaxation, 

TP'  °n°  obtalna  the  fast  absorber  equation  studied  by  Haus. 
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I I.  Some  Consequences  of  the  Basic  Equation 


We  have  not  found  a closed  form  solution  for  (5) . It  is 
possible,  however,  to  obtain  some  general  information  on  the  pulse 
shape  by  considering  the  significance  of  the  various  terms  in  (5) 
When  Ta  -*■  0,  it  was  shown  by  Ileus  that  6 + g = 0,  no  advance 
or  delay  of  the  pulse  takes  place.  This  is  clear  from  the  fact 
that  the  diffusion  (in  time)  of  the  pulse  as  caused  by  the  dis- 
persion preserves  the  symmetry  of  the  pulse,  and  so  does  an  absorber 
with  instantaneous  response.  When  f-  0,  it  is  possible  to 
determine  the  sign  of  6 + g by  the  following  reasoning:  an 

absorber  with  nonzero  relaxation  time  will  tend  to  "shave  off"  the 
front  part  of  the  pulse,  hence  causing  a net  delay,  6 + g < 0. 

In  the  wings  of  the  pulse  the  terms  proportional  to  v2 
may  be  neglected.  Hence  the  linear  equation  for  the  pulse  shape 
in  the  wings  is  given  by: 


1 + q 
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+ 6 * g A 

ul  dt 
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(10) 


The  solutions  of  this  equation  are  exponentials.  Hence  the  leading 
and  trailing  edges  of  the  mode-locked  pulses  are  exponentials,  exp(st) 
The  time  constants  of  the  exponentials  are  found  to  be: 
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In  the  front  end  of  the  pulse  the  growing  exponential 
is  the  solution,  the  tail  end  of  the  pulse  exhibits  the  decaying 
exponential.  Since  (6  + g)/2g  is  negative  and  1 + q - g > o 

for  stability  reasons,  one  finds  for  the  exponential  growth  factor 
s_  in  the  front  end. 


6 + g 


+ 1 + q - g - - -JL 


The  decay  factor  s+  < 0 of  the  tail  end  i: 


- 6 + g 


+ v(l TV  +1+q'9>  i«-i- 


V7e  find  that  the  pulses  must  be  asymmetric,  exhibiting  a slowly 
rising  leading  edge  and  more  rapidly  decaying  edge . 


A-ll 


III.  Computer  Solution 


Single  pulse,  periodic  solutions  to  Equation  (7)  were  found 
by  computer.  The  pulses  are  asymmetric  with  the  trailing  wing 
sharper  than  the  rising  wing  (see  Figure  2)  . Asymmetric  pulses 
of  this  typo  have  been  seen  experimentally  for  a system  in  which 
the  absorber  relaxation  time  was  shorter  than  the  inverse  gain 
linewidth'61.  The  pulsewidth  (full  width  at  half  intensity)  t„ 
is  found  to  be  only  a weak  function  of  absorber  relaxation  time" 
becoming  broader  with  increasing  T,,  as  shown  in  Fig.  3.  The 'eigen- 
value («  ♦ g)  as  a function  of  absorber  relaxation  time  is  shown 
in  Figure  4.  The  exponential  growth  and  decay  in  the  wings  is 
correctly  predicted  by  (12)  and  (13)  . 

The  absorber  loss  as  a function  of  time  during  the  passage  of 
a pulse  is  shown  in  Figure  -2  in  the  case  of  well-separated  pulses 
and  long  absorber  relaxation  time  (TA  » 10  tp>  . The  absorber 
loss  recovers  after  the  pulse  has  decayed.  However,  net  gain  exists 

« the  system  for  a critical  time  until  the  absorber  loss  relaxes 
back  sufficiently  to  cause  net  loss. 


m 
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IV.  Stability 


The  fact  that  there  is  net  gain  in  the  system  following  a 
pulse  in  the  case  of  long  nbsoibor  relaxation  time  may  lead  to 
instability  due  to  perturbation  growth.  The  evolution  of  a per- 
turbation introduced  into  the  steady  state  will  now  he  considered. 

The  rate  of  change  of  a pulse  over  many  passes  can  be  written 


— V(t,  T)  = 
8T 


- — 1 + q (t)  - q ( 1 + 

2 Q ' l 


3 j + 9 + 5 J)_ 


where  the  pulse  count  n has  been  replaced  by  a continuous  time 
variable,  T = n-TR.  The  two  time  variables  are  assumed  indepen- 
dent of  each  other;  T measures  the  pulse  repetitions,  and  t 

measures  the  "local"  tine  within  a pulse.  In  the  steady  state, 

( 3/9T) v ( t , T)  = 0. 

To  investigate  stability  we  introduce  a perturbation  6v  to 

the  steady  state  pulse  solution  vQ.  The  waveform  in  a period 

specified  by  T becomes  (compare  Ilaus  • analysis  in  Section  IV  of 
Ref.  (3)) 


v(t,  T)  = vQ(t)  + 6 v ( t , T)  . 


The  perturbation  is  allowed  to  change  on  the  slow  time  scale  T. 
In  the  region  where  the  perturbation  does  not  overlap  the  steady 
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•slate  solution  its  influence  on  the  absorption  and  gain  can  be  neglected. 
V.hcre  overlap  occurs  the  added  saturation,  must  be  included  if  the  per- 
turbation is  in  phase  with  the  steady  state.  Stable  operating  regimes 

with  respect  to  inphase  perturbations  have  been  found  by  Haus  for  the 
fast  absorber. 

In  the  case  of  tile  fast  absorber,  inphase  perturbations  of  a stable 
operation  decay,  because  the  decrease  of  absorption  caused  by  a perturba- 
tion is  compensated  by  a larger  decrease  of  gain.  Quadrature  perturba- 
tions are  neither  stabilized  nor  destabilized.  It  is  reasonable  to 
assure  that,  in  the  present  case  of  a slow  absorber,  the  steady  state 
pulse  is  stable  with  respect  to  inrhasc  perturbations  overlapping  with 
the  pulse,  on  the  other  hand,  it  is  clear  that  all  perturbations  fol- 
lowing the  steady  state  pulse  and  not  overlapping  with  it  grow  because 
they  do  not  affect  the  gain  as  produced  by  the  steady  state  which 
is  positive  for  some  tire  following  the  pulse  (wo  assume  the  saturable 
absorber  relaxation  time  is  long  enough  that  the  net  gain  fol- 

lowing the  pulse  persists  for  at  least  several  pulsewidths)  . As  the 
perturbation  grows,  it  is  pulled  progressively  into  overlap  with  the 
steady  state  pulse.  Perturbations  in  quadrature  with  the  steady  state 
pulse  will  be  least  stabilized  because  they  affect  neither  the  gain 
nor  the  loss  when  they  core  into  overlap  with  the  steady  state  pulse. 

Our  objective  is  to  determine  the  maximum  growth  possible  for  quadra- 
turc  perturbations . 

The  evolution  of  the  perturbation  in  T is  governed  by  the  lo.ss 
modulation  g(t)  established  by  the  steady  state  pulse.  In  general, 

the  perturbation  experiences  three  major  changes  on  the  slow  time 
scale : 


(1)  Amplitude  growth  due  to  the  excess  gain, 


A-14 


(2)  motion  toward  the  steady  state  pulse  due  to  the  fact  that 
the  steady  state  pulse  experiences  a delay  in  each  round  trip,  and 

(3)  width  variation  due  to  the  gain  dispersion  and  loss  modula- 
tion curvature. 

Ultimately  the  perturbation  moves  into  the  steady  state  pulse. 

If  we  assume  a gaussian  initial  perturbation 


6v  = a^  exp 


(t  - t.)2 
~ 2t  . 2 

l 


(16) 


a closed  form  solution  to  (12)  for  the  evolution  of  the  perturbation 
is  possible  provided  we  divide  the  loss  modulation  following  the 
steady  state  pulse  into  two  regions  as  shown  in  Figure  5.  The  ex- 
cess gain  is  approximated  by  a straight  line  in  Region  I,  where  the 
loss  is  relaxing  back  after  being  depleted  by  the  pulse;  and  by  a 
parabola  in  Region  II,  where  the  loss  is  being  depleted  by  the  steady 

state  pulse.  The  results  of  the  analysis  are  presented  below.  De- 
tails appear  in  Appendix  A. 

We  neglect  the  change  in  perturbation  width  during  its  evolution. 
For  a perturbation  having  an  initial  width  approximately  equal  to  that 

of  the  steady  state  pulse  and  much  greater  than  the  inverse  laser  line- 
width*,  i ,e  . 


T . - t >> 
i s 


2w, 


(17) 


the  expression  for  the  total  amnl  i fade  growth  of  the  ;cn 

initiated  at  t^  = T^/2  assumes  the  form 


* thi5  i£:.not  contradictory  since  the  steady  state  pulst 

ur.es  only  a fraction  of  the  laser  linowidth. 
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a = a^  exp 
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The  expression  is  separated  into  factors  which  describe  the  growth 
m the  two  separate  regions  of  the  loss  modulation  following  the 
steady  state  pulse.  As  might  be  expected  the  growth  in  region  I is 
directly  dependent  on  the  absorber  relaxation  time  , while  that 
m region  II  is  only  dependent  on  via  the  steady  state  "delay" 

(g  + 6)  and  pulsewidth  ts  which  are  plotted  in  Figures  4 and  5. 

The  approximate  dependence  of  perturbation  growth  on  absorber 
relaxation  time  given  by  equation  (16)  is  illustrated  in  Figure  6. 
Also  plotted  is  a computer  solution  for  a gaussian  perturbation 
having  an  initial  width  equal  to  that  of  the  steady  state  pulse. 

The  difference  between  the  computer  and  approximate  curves  for 
large  results  from  our  neglect  of  shaping  during  the  growth 

of  the  perturbation  in  the  derivation  of  (16)  . 

If,  for  example,  we  set  as  a criterion  for  stability  that  the 
perturbation  amplitude  at  no  time  exceed  one  tenth  that  of  the  steady 
state  pulse,  then  for  an  absorber  relaxation  time  ytj  . 5> 
equation  (16)  requires  the  initial  perturbation  amplitude  a,  to 
be  at  least  40  du  down  from  the  steady  state  pulse  amplitude. 
Fundamental  noise  in  the  system,  o.„.  spontaneous  emission  noise. 
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I 

certainly  satisfies  this  constraint.  Other  sources  of  noise,  such 
as  scattering  from  intracavity  apertures,  might  lead  to  instability. 

Thus,  stable  steady  state  nodelocking  with  slow  saturable  ab- 
sorbers is  possible  for  some  regime  of  > t , without  action  of 
the  laser  medium  as  required  for  a truly  slow  saturable  absorber ^ 

>:>  Tp^  * Stability  is  made  possible  by  the  fact  that  perturbations 
appearing  in  the  system  are  not  free  to  grow  without  bound,  but  are 
constrained  by  the  loss  modulation  of  the  absorber  to  approach  a 
steady  state  amplitude. 
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V . Conclusion 


Vto  ha vo  obtained  steady  state  single  pulse  solutions  for  the 
problem  of  mode locking  with  a slow  saturable  absorber  where  the 
gam  in  assumed  constant  over  a pulse  repetition  time  and  the 
pulsewidth  is  much  less  than  a pulse  repetition  time.  The  pulses 
were  found  to  bo  asymmetric  with  a slowly  rising  edge  and  more 
rapidly  decaying  edge.  The  modelocked  pulsewidth  is  a weak  func- 
tion of  the  absorber^ relaxation  time.  In  contrast  to  the  assertions 
of  previous  authors151  v,e  have  found  that  stable  modelocking  is 
possible,  even  in  cases  where  the  absorber  relaxation  time  Tft  is 
comparable  to  or  greater  than  the  pulsewidth.  since  the  steady 
state  pulsewidth  utilises  only  a fraction  of  the  laser  linewidth, 

stable  operation  Is  possible  even  for  u T >>  l. 

LA  * 
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Appendix  A 


of  a Perturbation  Pol  loving  the 


Steady  State  Pul se 


AS  illustrated  in  figure  5 the  loss  modulation  following  the 
Steady  state  pulse  is  divided  into  two  regions  in  whieh  the  excess 
gam  is  approximated  as  follows: 


Agj  - Agm  ^1  - JL  j 


4,I1  = l1  ' — 


t ^ 

II 


0 1 fc  1 tI 


~tII  1 t < 0 . 


(A  .1 ) 


(A. 2) 


he  paramo. cis  Agm,  tj  and  t^  characterize  the  magnitude  and 

exceos  gam  respectively  and  are  approximately  related 
to  the  steady  state  system  parameters  as  follows: 


Ag  = 1 + 


q - g 


(A. 3) 


tl- !!s 

1+q-g+Ag  A 2 

■ m 


(A  .4 ) 


t„  = 


4 In 2 V g Wl  Ts  * 


(A. 5) 


Introducing  the  spore* i motion  fn  *-?>«  *•  ■: ™ 

‘ - ..  m rr,  t;„m  time  dependence  of  the  ex  cos: 


Win  (A.  1)  and  (A. 2)  into  the  equation  of 


motion  obeyed  by  the 
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perturbation  (10)  , wc  obtain 


— <Sv(t,  T)  = 
9T 


20 


Ag(t)  + _S il  + I g + A | j_ 


“l*  9t2 


WL  9t 


ov.  (A. 6) 


We  further  assume  the  perturbation 


to  have  a gaussian  shape 


<5v(t,  T)  = ao(T)  exp 


[t  - tQ  (T)  ] 2 


2tq2 (T) 


(A. 7) 


characterised  by  a nlowly  varying  aKputudo  . 
and  width  to(t).  Initial  condition 

to  be 


0(T),  timing  t (T)  . 

o 


on  the  perturbation  are  taken 


ao(°)  =.a. 


V°>  = t. 


V(0»  ” V 


ThC  °f  <A-l)  ami  (A'7)  in  <*■*>  load,  to  the  following  act  o, 
equations  governing  the  evolution  of 

amplitude  in  Region  I. 


perturbation  width,  timing  and 


dx  w 

O _ _^0  g 


dT  20  wt2t 
L o 


(A.  9) 


! 
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dt  to  | 

o = | g + 6 

dT  2Q  Wr 


Aym  to- 
rn L 

:i  ' | g + 6 1 ° 


(A. 9) 


— = ~ Ag  (l  - 2 

dT  2Q  [ \ t '/  ,,  *T  2 

hO 


(A. 10) 


If  wo  restrict  our  attention  to  cases  where 


—2—  -2.  T « 1 
wl  2 T i 2 Q 


(A. 11) 


in  Region  I,  then  the  width  variation  of  the  perturbation  can  be 
neglected  and  solutions  for  timing  and  amplitude  are  found  to  be 


to(T)  . t.  - ^ + T 

fcI  \ AVi2  WL  1 Q 


(A. 12) 


ao  (T) 


i eXP  A^m  1 - 


t.  \ w 
1 | o 


t i/2Q 


g + 6 


Ag  t.  - 
m i 


(A. 13) 


The  perturbation  reaches  the  boundary  of  Region  I,  t =0 

' o ' 

at  time  T = t given  by 
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^ T * [ ]_  + J_q  + 6 1 


Agt  t . 2 
i 


(A. 14) 


Thus,  a perturbation  initiated  a<-  f _ < . 

-nil rated  at  tf  - tj.  has  an  amplitude  at 


4 = 0 


a^  exp 


l0L  4 


ip 

1 + !il 

tj  | g + 6 


(A. 15) 


Th°  PCrtU‘1,ntion  enter=  Region  II  witl,  the  initial  con- 
ditions to  . t.,  ao  . alt  Subsequent  evolution  of  the  perturbation 
is  governed  by  equation  (A.6)  where  (A.2)  is  used  to  approbate  the 

excess  gain.  To  simplify  the  analysis  we  shall  assume  that  the 

perturbation  width  on  enterina  Roninn  tt 

ring  Jtegion  II  is  approximately  equal  to 

the  steady  state  width,  i.  e. 


Ti  * Ts 


(A. 16) 


SO  that  changes  in  width  can  be  neglected  throughout.  The  equations 
governing  timing  and  amplitude  evolution  then  become 


0) . | . / \ 
- - -a.  ls-i-il  (i  + i ) 

20  “l  V 6S  / 


(A. 17) 


\ 
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da 


£ = “o  | g + <S|2 
dT  2Q  4 g 


t 


1 - 


(A.  10) 


where  the  steady  state  timing  approached  by  the  perturbation  i: 
that  of  the  steady  state  pulse,  namely 


ts  = - -L5L+  M / 4l  X2 
2 WL  \Ti 


(A. 19) 


Equation  (A. 17)  is  of  the  form  of 
tion  rate 


a relaxation  equation  with  relaxa- 


a - Ag, 


m 


^t/  q 


(A. 20) 


Hie  slow  time  evolution  of  t i 


D is  given  by 


t 

o 


“ V1 


(A. 21) 


Integrating  (A. 18)  then  gives 


a(T) 


aj  exp 


4 /Kq~  g ^ L 11 


(A. 22) 


in  the  limit  where  T >>  1/ct. 


V?i-h  the  use  of  ;A.3)-(«.5)  equations 


t * m 

C**  i va 


/ » on 


eujuD.tnc 


to  give  equation  (14)  of  the  text.  Assumptions  (A. 11)  and  (A.16) 
restrict  the  validity  of  the  analysis  presented  hero  to  cases  where 
the  perturbation  width  variation  can  be  neglected.  The  inclusion 
of  width  variation,  however,  loads  to  a smaller  growth  than  predicted 

by  (14),  so  that  (14)  represents  an  upper  bound  on  the  possible 
perturbation  growth  . 


L_ 
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Figure  Captions 


•—?ure JL*  Passively  Fodelocked  Laser  System 


Absorber  Loss  as  a Function  of  Time  for  a Well- 
Separated  Pulse  with  Long  Absorber  Relaxation  Time 
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Transient  Buildup  of  Fast  Saturable  Absorber  Modelockin 

by 


Christopher  P.  Ausschnitt 
Massachusetts  Institute  of  Technology 
Research  Laboratory  of  Electronics 

Abstract 

In  the  following  paper  we  present  a closed  form  deterministic 
solution  to  the  buildup  of  fast  saturable  absorber  modelocking  in 
the  homogeneously  broadened  laser.  The  evolution  of  the  modelocked 
waveform  from  one  transit  through  the  laser  cavity  to  the  next  is 
treated  adiabatically,  as  a succession  of  guasisteady  states.  The 
secant  hyperbolic  pulse  solution  found  by  Haus  is  reached  in  the 
limit  of  widely  separated  pulses. 
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Introduction 


Previous  authors  have  taken  a statistical  approach  to  the 
study  of  the  buildup  of  passive  modelocking, f1"51  stating  that 
the  occurrence  of  modelocking  relies  on  the  selective  bleaching 
of  the  absorber  by  the  highest  noise  spike  in  a cavity  transit 
time.  Although  it  is  true  that  lasing  initiates  with  spontaneous 
emission  noise,  two  facts  argue  against  the  statistical  approach: 

1*  During  the  initial  buildup  of  radiation  in  the  cavity, 
where  saturation  of  the  absorber  Js  negligible,  the  bandwidth  of 
the  noise  is  narrowed  in  successive  passes  through  the  cavity 
due  to  the  gain  dispersion. 

2.  Passive  modelocking  is  consistently  reproducible  if  the 
laser  is  operated  near  threshold  and  careful  control  of  the  system 
parameters  is  maintained.^ 

The  first  of  the  above  observations  indicates  that  the  initial 
spontaneous  emission  noise  is  smoothed  by  gain-narrowing  so  that 
the  sharp  noise  spikes  are  "washed  out".  In  fact,  the  near  thres- 
hold operation,  stated  in  the  second  observation,  limits  the  in- 
itial bandwidth  of  the  gain  medium,  i.e.  the  number  of  cavity  modes 
which  see  net  gain.  Thus,  near  threshold  operation  ensures  that 
gain  narrowing  will  "smooth  out"  the  noise  fluctuations.  We  argue, 
therefore,  that  by  the  time  the  absorber  saturation  becomes  signi- 
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ficant  the  oscillation  state  in  the  laser  can  be  described  by 
a perturbed  single-cavity  mode  (SM)  quasi-steady  state. 
Consequently,  the  buildup  of  modelocking  cannot  be  due  to  the 
selective  bleaching  of  the  absorber  by  a noise  spike. 

We  show  here  that  modelocking  results  from  the  fact  that  at 
a definite  power  level  in  the  cavity  the  SM  state  to  which  the 
initial  noise  evolves  becomes  unstable  with  respect  to  a sinusoidal 
perturbation.  Thus,  an  expression  for  the  modelocking  threshold, 
defined  as  the  point  at  which  the  SM  acquires  a sinusoidal  modu- 
lation is  found.  The  subsequent  evolution  of  the  system  from  SM 
to  steady  state  modelocked  operation  is  deterministic  and  is  treated 
as  a succession  of  the  steady  state  solutions  to  the  modelocking 
equacion  derived  by  Haus  . [ ?]  Furthermore,  for  each  set  of  system 
parameters— characterizing  the  cavity,  saturable  absorber  and  gain 
medium— a unique  stable  periodic  pulse  solution  is  reached.  In 
the  limit  of  well  separated  pulses  the  steady  state  solution  is 
that  of  the  single  secant  hyperbolic  pulse  found  by  Haus. 


A-35 


— — ?he  Single  Pulse  Steady  State  Modelocking  Solution 

If  the  mode  locked  laser  system  consists  of  a short  section 
of  saturable  absorber  and  a short  laser  gain  medium,  and  under  the 
assumption  of  weak  absorber  saturation,  the  equation  governing  mode- 
locking  by  a fast  saturable  absorber  has  been  shown  by  Haus  to  be:17 

_ 2 8v(t,  n) 

4“c  tr  Sn 

(1.1) 


1 + q - g 


_ q| v(t,  n) | 2 
PA 


to. 


3t‘ 


v(t,  n) 


where 


v = field  envelope 

t = time  "local"  to  a single  cavity  round  trip  transit  time 

n = integer  which  counts  the  number  of  round  trip  transits, 

treated  as  a continuum  variable  over  many  transits 

tr  = round  trip  transit  time  in  cavity 
Awc  = cavity  ^ode  linewidth 

q = small  signal  absorber  loss 

PA  = saturation  power  of  the  absorber 
g = saturated  gain 
WL  = laser  medium  half  width 


The  saturated  gain  is  related  to  the 


power  in  a single  round  trip 


I 
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transit  by 


g(n) 


Vn> 

P(n) 


(1.2) 


where  gQ(n)  is  the  known  small  signal  gain,  PL  is  the  saturation 
power  of  the  homogeneously  broadened  laser  medium  and  the  power  is 


P (n)  = — 


V2 


tr  J -v2 


v(t , n) 1 2 dt. 


(1.3) 


It  should  be  noted  that,  in  contrast  to  the  fast  absorber  equation 
used  by  Haus,  we  have  neglected  the  linewidth  broadening  due  to 
gain  saturation  by  making  the  substitution 


-2—  = L.t  <3  = JL 


The  approximation  is  justified  by  the  fact  that  the  analysis  is  only 
valid  in  the  regime  where  the  gain  is  near  threshold.  Thus,  gain 
variation  is  negligible  compared  to  the  normalized  cavity  loss, 
but  is  significant  compared  to  1 + q - g.  Although  the  following 
analysis  can  be  carried  out  without  making  this  approximation,  the 
inclusion  of  power  broadening  sacrifices  simplicity  for  little 
additional  insight  into  the  problem. 


.j 
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Also,  many  systems  are  bandwidth  limited  by  a dispersive 
element  (e.g.  etalon)  inside  the  laser  cavity  in  which  case  the 
diffusion"  operator  expressing  pulse  broadening  in  time  is  gain 
independent. 

In  the  steady  state  the  waveform  v experiences  no  change 
from  one  transit  through  the  cavity  to  the  next.  Thus,  8v/3n  = 0 
and  equation  (1.1)  reduces  to 


1 + q - g - 


_1 

— L? 


v(t)  = 0. 


(1.4) 


Haus  has  shown  that  one  stable  solution  to  (1.4)  is  a single 


pulse 


v(t)  = 


(1.5) 


cosh  — 


where  the  pulse  parameters  vq  and  x are  determined  by  the 

r 

system  parameters  via  the  eigenvalue  relations 


1 + q - g = SL  = _ 1 
2 


(1.6) 


in  conjunction  with  (1.2).  In  the  following  section  we  show  that 

the  general  solution  to  (1.4)  is  periodic,  and  that  (1.5)  holds 

only  in  the  special  case  where  T >>  x . 

R p 
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— — — - PGriodic  Steady  State  Mode  locking  Solution 

The  secant  hyperbolic  steady  state  solution  found  by  Haus 
is  only  one  solution  to  equation  (1.4).  The  fact  that  it  is  a 
single  pulse  or  "soliton"  solution  implies  that  it  corresponds 
to  the  case  where  the  cavity  transit  time  is  much  longer  than  the 
pulsewidth , TR  » Tp,  so  that  the  inherent  periodicity  of  the 
system  can  be  ignored.  During  the  buildup  of  modelocking  from 

SM/  however'  ifc  is  dear  that  the  periodicity  of  the  system  must 
be  taken  into  account  since  the  pulsewidth  Tp  is  initially  com- 
parable  to  TR.  As  a first  step  toward  understanding  the  transient 
buildup,  therefore,  we  must  find  the  periodic  solution  to  (1.4). 

In  general,  equation  (1.4)  is  the  equation  of  motion  of  a 
particle  in  the  potential  well 


U(v)  = 


-3Xl+l 

4P„  2 


(1  + q - g)  v2  + C 


(2.1) 


where  C is  a constant.  If  the  particle  is  launched  with  zero 
velocity  at  the  well  height  0,  corresponding  to  a displacement 
V°'  osclHates  between  the  turning  points  defined  by  the  roots 
of  U(v).  Thus,  the  solution  is  periodic,  with  time  dependence 


v(t)  = v ll  - y sn : 


(r 


(2.2) 
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where  sn  is  a Jacobian  elliptic  function  of  t/t  of  modulus 

[ 8 ] P 
Y and,  for  convenience,  we  have  defined  the  new  constant  Y as 


q v ** 

Y = — C + 1. 

4PA 


(2.3) 


Substitution  of  (2.2)  in  (1.4)  gives  the  eigenvalue  relations 


2PA 


(2.4) 


and 


1 + q - g = — (2  - y ) (2.5) 

' it  Tp 

which  determine  the  peak  amplitude  v and  "pulsewidth"  t 

° p 

within  each  period.  The  periodicity  of  the  solution,  which  we 
require  to  equal  the  cavity  transit  time  TD,  determines  the 
constant  y via  the  relation 

2K(*>  ■ ^ (2.6) 

where  K is  the  complete  elliptic  integral  of  the  first  kind  of 
argument  y.  In  the  limit  TR  >>  xp,  equation  (2.6)  dictates 
that  Y - 1 and  equations  (2.2  - 2.4)  become  the  secant 
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hyperbolic  solution  and  eigenvalue  relations  found  by  Haus  (1  5 
’*  ReqUiri"9  ^ TOOtS  to  be  real  sets  the  limits 

on  Y/ 


0 < Y < 1 


(2.7) 


.. ...» .... 

Tp  via  (2.6)  , namely  R 


_R 

P 


Tr  < < co, 

T 


(2.8) 


The  lower  limit  nf  /?  o\ 

fora,  v ■ corresponds  to  the  case  where  the  wave- 

, 1S  SH  “ thtinitesimal  sinusoidal  ripple  one  hil 

wavelength  of  which  is  fitted  in 

earlier  th  ‘ tOU"d  trip  tra"-t.  As  noted 

earlier  the  upper  limit-  ; „ .. 

P "it  is  the  case  where  v is  a single  se  -ant 

hyperbolic  pulse.  nt 

The  power  in  a sinale  franco 

tutina  th  “ determined  ^ substi- 

tuting the  general  solution  (2  2)  in  „ 

of  |v„  ' -2>  equation  (1.3).  The  integr, 

I over  one  period  is  expressible  in  tents  of  E,  the 
complete  elliptic  inteorai 

„„  9 °f  the  second  kind  of  argument  v 

so  that  (1.3)  gives  Y' 


2vqJ  Tp 

T« 


e(y)  = V 2 Sill 
° K(y) 


(2.9) 


r - 
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III.  The  Modelocking  Threshold 

Lasing  initiates  when  the  small  signal  gain  gQ  first  breaks 
above  the  small  signal  loss  1 + q.  In  the  region  near  the  onset 
of  lasing  we  assume  that  the  saturation  of  the  absorber  is  negli- 
gible. The  initial  bui 3 of  the  field  is  described  by 


_2 3v(n,  t) 


A“cTR 


" 9 -(l  + <l)  + -1-  -4  V(n,  t). 

L &L  3t J 


(3.1) 


We  can  express  the  field  in  the  nth  transit  by  a superposition 
of  the  cavity  modes,  namely 


v(n,  t)  * £ Vk  exp  j + /"  Bk  (n)  dn 

. L TR  J 


(3.2) 


where  k is  an  integer  which  labels  tne  modes,  VR  is  the  am- 
plitude and  3k  is  the  gain  coefficient  characteristic  of  the 
kth  mode- 

Substitution  of  (3.2)  in  (3.1)  gives  an  expression  for  the 
gain  coefficient 


AwcTR 


Bk(n)  - g(n)  - (1  + q)  - -L-  (k  JL.V  . 


(3.3) 


The  mode  nearest  line  center  (k  - 0)  experiences  the  most  rapid 
growth  due  to  the  band  narrowing  effect  of  the  gain  dispersion. 
Thus,  operation  on  the  single  mode  (SM)  at  line  center  win 


at  line  center  will 


A-42 


dominate.  We  assume,  therefore,  that  the  laser  prior  to  the  onset 
of  modelocking  is  m EM  operation,  namely  that  its  output  does 


not  vary  over  the  transit  time  T 

R 


Thus  far  we  have  neglected  the  saturation  of  the  absorber. 


Ultimately,  however,  the  SM  described  by  v (n)  grows  to  the 


point  where  the  equations  describing  the  subsequent  growth  of  the 
field  must  include  the  absorber  saturation.  Hence,  the  equation 
governing  growth  becomes 


dv 


dn 


" Mn>  v (n) 


(3.4) 


where  8Q  is  given  by 


2Bo(n)  , , q v 2 (n) 

■ = g(n)  - (1  + q)  + 2 


TR 


(3.5) 


A 


Note  that  gam  dispersion  plays  no  role  in  determining  & the 
gain  coefficient  of  the  mode  at  line  center.  Since  the  field  is 
"time  independent"  within  T the  saturated  gain  is 


q (n)  = 


go  (n) 


1 + 


(n) 


(3.6) 


The  modelocking  threshold  is  determined  by  the  value  of 


Vo (n)  for  which  the  first  modolocking  solution  to  (1.4)  is  ob- 
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tained.  As  shown  by  (2.8),  tho  first  solution  is  the  case  where 
the  EM  waveform  acquires  a sinusoidal  ripple,  with  one  half 
wavelength  fitted  in  a round  trip  transit.  Thus,  the  modelocking 
threshold  is  defined  by  n = £ such  that 


R 


TpU) 


(3.7) 


and  we  have 


v*>  - 


u. 


,TR  V q 


(3.8) 


The  SM  field  must  build  up  to  a value  ,„(*)  before  blocking 
initiates . 

To  prove  that  the  transition  from  the  SM  to  the  first  mode- 
locking solution  is  deterministic,  we  now  show  that  a sinusoidal 
perturbation  on  the  SM  experiences  more  rapid  growth  than  the 

SM.  in  other  words,  the  system  must  modelock  rather  than  continue 
to  run  in  a single  mode. 

Returning  to  equation  (1.1),  we  assume  a perturbed  solution 
of  the  form 


v(n,  t)  - v (n)  + r {v  a 

k 5/  0 k 


jk  + 
- tr 


3^ (n) dn 


(3.9) 
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'■■>eve  vQ<n)  is  the  growing  SM  solution  whose  gain  coefficient 
8o  9^vsn  by  (3.5) . The  perturbation  has  been  expanded  in  the 
cavity  -nodes.  Substituting  (3.9)  in  (1.1)  „e  find  that  to  first 


order  0k(n)  is  given  by 


Aw  T_ 
c R 


Bk(n)  = g(n)  - (1  + q)  + i£L  v 2 (n) L.  ( k 2jt 


(3.10) 

Since  the  perturbation  is  "orthogonal”  to  the  single  mode  field 
over  a cavity  transit  time,  the  saturated  gain  does  not  change  to 
first  order  in  the  perturbation  amplitude,  and  is  still  specified 
by  (3.6) . From  (3.5)  we  see  that  the  difference  between  the  gain 
coefficient  of  the  single  mode  field  and  that  of  the  perturbation 


(Bk  - 0 ) «=  13.  v 2 

Aw  T K 0 P o 


k 21 


(3.11) 


Maximum  growth  is  experienced  by  the  perturbation  for  which  k = 1. 
Setting  Bj  = then  extablishes  the  value  of  vQ  for  which  the 
perturbation  growth  first  exceeds  the  SM  growth.  Clearly  the 
threshold  value  of  vq  predicted  by  (3.11)  is  the  same  as  that 
given  by  (3.8).  Thus,  for  vo>vo(f)  the  sinusoidal  perturbation 


k 
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grov/s  faster  than  the  SM.  Consequently,  once  the  threshold 

vo  c voU)  reached,  the  system  must  ease  into  the  modelocking 

solution . 


J 


IVj — The_Transient  Buildun 
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We  have  shown  in  the  previous  section  that  at  a particular 
SM  freld  amplitude  in  the  cavity  the  single  mode  field  must  break 
into  oscillation,  furthermore,  the  initial  oscillation  corresponds 
to  the  first  steady  state  solution  of  the  modelocking  equation. 

The  transient  evolution  from  the  initial  oscillation  to  the  final 
steady  state  ir  now  treated  as  succession  of  quasi-steady  state 
solutions  to  equation  (1.1),  where  the  left  hand  side  of  the  equa- 
tion dictates  a change  in  the  waveform  v(„,  t,  from  one  transit 
to  the  next. 

To  determine  precisely  how  the  waveform  evolves  we  constrain 

the  buildup  to  be  adiabatic,  namely  the  increase  in  the  field  energy 

transit  is  set  equal  to  the  energy  supplied  by  the  system  in 

one  transit  so  that  the  total  energy  remains  constant.  The  energy 
balance  equation 


4Wfield  = - AW 


system 


(4.1) 


oan  be  obtained  directly  from  (1.4)  by  multiplying  both  sides  by 

v*(n,  t)  and  integrating  over  one  period  or  transit  tine.  Thus, 
(4.1)  becomes 

— y M!  dt  - aocTR  Ug  - (1  + q)  )/ I v|  s dt  + JL  /|v|*  dt 

D J 


1 f dvl 2 . 

7 r "H  dt} 

fit  1 J ^4. 


WL  • dt 


(4  .2) 
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where  f denotes  integration  over  one  period  - Tr/2  < t < t /2  . 

Ke  assume  that  the  power  in  the  field  experiences  a charge 
over  one  transit  time  given  by 


P (n  +!)  = [!  + 2a(n)]  P(n) 


(4.3) 


where  a(n)  is  the  growth  rate  of  the  field  to  be  determined 
from  the  energy  balance  equation.  Thus,  for  two  energy  states 

which  are  differentially  close,  the  left  hand  side  of  the  energy 
balance  equation  becomes  simply 


cl  f 

AWfield  = “^Jiv|2dt  = 2tr  Q(n)  P(n)  . 


(4.4) 


m order  to  evaluate  the  change  in  the  system  energy  we  make 

use  of  the  fact  that  the  quasi-steady  state  solution  for  each  value 
Of  n is  given  by  (2.2),  namely 


v2  (n,  t)  = ini ) 
E (y  in)  ] 


p(n)  1 - y (n) 


Tp(n) 


Y (n) 


(4.5) 


Tp2  (n)  g j.  P^ 

WL2  9 PL  K fY  (n)  ] P(n) 


(4  .6) 
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i 


** 


* 


The  constant  y(n)  and  the  saturated  gain  g(n)  are  determined 
by  the  relations 


Tr 

2K[y(n)  ) = — 

Tp(n) 


and 


g(n) 


Vn> 

i + ' 


(4.7) 


(4.8) 


Substitution  of  (4.5)  in  the  right  hand  side  of  (4.2)  and  equating 
to  (4.4)  gives  the  desired  expression  for  a(n) 


— - — ft  (n)  = [g (n)  - (1  + q)  J + [2  - y(n))  — iiY.fo.IJL  # 

AucTr  2E  [y  (n)  J PA 

(4.9) 

The  equations  governing  the  transient  buildup  are  now  com- 
plete. The  waveform  v(n,  t)  for  each  value  of  n is  given  by 

(4.5)  where  the  pulsewidth  in  each  period  t (n)  , the  constant 

P 

Y (n)  and  the  saturated  gain  g(n)  are  determined  by  equations 
(4.6  - 4 . 8)  . The  evolution  of  the  nth  waveform  to  the  (n  + l)tn 
waveform  is  dictated  by  equation  (4.3),  where  conservation  of 
energy  constrains  q(n)  to  be  given  by  (4.9). 
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The  initial  condition  on  P(n)  was  shown  in  section  III 


to  be 


_ 2 2P 

P (in)  = -I t 

wrTn  q 


(4.10) 


The  final  steady  state  defined  by 


v(t,  n + 1)  = v { t , n) 


(4  .11) 


requires  that 


fff  = 5<n)  - (1  + + [2  - Y {n , J . 0. 

CTR  2PAK[y(n)] 


(4.12) 


We  recognize  that  (4.12)  is  the  eigenvalue  relation  (2.5)  determined 
earlier  to  describe  the  general  steady  state  solution.  Equations 
(4.12)  and  (4.6  - 4.8)  now  uniquely  specify  the  final  steady  state 
solution  as  was  shown  in  section  II . 

Of  course,  m cases  where  the  small  signal  gain  of  the  laser 
medium  never  reaches  a constant  value  (e.g.  gain  or  Q switched 
lasers)  condition  (4.11)  is  never  satisfied.  This  does  not  limit 
the  effectiveness  of  the  analysis  provided  the  gain  variation  over 
the  cavity  transit  time  TR  is  small.  In  fact,  the  analysis  en- 
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ables  us  to  trace  the  complete  evolution  of  transient  modelocking 
from  start  to  finish. 
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